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Abstract — This paper proposes a simple approach to design 
a PID controller for an uncertain SISO system that will 
guarantee closed loop stability when the input is varied 
between and 60 % of its working range. The Robust 
Controller design method used here is based on Small Gain 
theorem. A family of transfer functions corresponding to 
different working ranges is used to describe the uncertain 
plant for which the controller is to be designed. In this paper, 
the controller is designed for a warm air-drying chamber 
which is a good example for an uncertain system. The design 
is carried out considering the plant behaviour at 3 working 
points. The closed loop step responses obtained, prove the 
robust operation of the designed controller for the plant 
under the considered working range. 

Keywords — Robust controller, Uncertain model, Small 
gain theorem, Closed loop stability. 

I. Introduction 

Any physical system can never be described by an exact 
mathematical model. The mathematical model is always only 
an approximation of the true system dynamics since we ignore 
certain factors such as non-linearities, time delay, effects of 
reduced order model, variations in system parameters due to 
environmental changes etc. which may affect the plant 
behaviour. This results in modelling errors and hence improper 
design of a controller. 

The deviation of the actual plant model from the nominal one 
is called model uncertainty. Uncertainty can be classified as 
disturbance signals and dynamic perturbations. The former 
includes input and output disturbances and the latter represents 
discrepancy between the mathematical model and the actual 
dynamics of the system. 

Unstructured uncertainty: 

Dynamic perturbations such as due to unmodelled high 
frequency dynamics can be lumped together into a single 
perturbation block A. This is called unstructured uncertainty. 
The unstructured dynamics uncertainty can be represented in 
different ways as shown in figures 1 and 2. G p (s) denotes the 
actual perturbed system dynamics and G (s) denotes the 
nominal model of the system. The block A is uncertain but is 
normally bounded. The additive uncertainty representations 
give an account of the absolute error between the actual 
dynamics and the nominal model. The multiplicative 
uncertainty representations indicate the relative errors. 



G p (s) = G (s) + A(s) 



Fig. 1 Additive Perturbation Model 



G p (s)= G (s)[l + A(s)} 



-M+ 



Fig. 2 Multiplicative Perturbation Model 
Structured uncertainty: 

An uncertain model can be described mathematically by a 
transfer function having uncertain parameters as its 
coefficients. These coefficients vary about their nominal 
values within the uncertainty bounds under different working 
conditions. 

Parametric uncertainty: 

In many systems dynamic perturbations may also be caused by 
inaccurate description of component characteristics. Such 
perturbations may be represented by variations of certain 
system parameters over some ranges. These affect the low 
frequency range performance and are called parametric 
uncertainties [i]. 

PID controllers are widely popular in industries because they 
are easy to implement. However they lack generality since 
they use heuristic tuning methods like Ziegler-Nichols [ii]. 
Dubravska M. and Harsanyi L. [iii, iv] have proposed a novel 
and easy to use method of designing a PI/PID controller for 
uncertain systems using Small gain theorem. Ying J. Huang 
and Yuan-Jay Wang [v] have demonstrated the design of a 
robust PID controller for an uncertain plant with bounded 
parameters using Kharitonov theorem. Clarke et al. [vi] 
introduced the Generalized predictive control to overcome the 
limitations of Model predictive control like open loop system 
being unstable, unknown plant orders etc. The computational 
burden has been significantly reduced by the GPC technique 
for FOPDT model proposed by Bordons and Camacho [vii, 
viii] . A control problem is formulated as a mathematical 
optimization problem in H M technique [ix, x, xi]. This paper 



IJSET@2014 



Page 207 



International Journal of Scientific Engineering and Technology 
Volume No.3 Issue No.2, pp : 207 - 210 



(ISSN : 2277-1581) 
1 Feb 2014 



USE! 

aims to use a simple, frequency domain based approach to 
design a PID controller for a third order uncertain system that 
is described by a family of Nyquist plots. 

II. Methodology 

Consider a SISO uncertain system G(s) described by a family 
of transfer functions and being controlled by a controller C(s). 
Then the loop transfer function is given by L(s) — C(s) G(s). 



R(s) 







U(s) 






C(s) 




G(s) 




► 



Y(s) 



Fig. 3 Closed Loop Control System 

According to Small Gain theorem, if the magnitude of the loop 
gain is less than one then, the closed loop system is robustly 
stable. Hence for closed loop robust stability, | 
L(s) | < 1 for all co €[0, oo] 

Suppose that the plant to be controlled is uncertain and 
described by a number of stable transfer functions. Let its 
nominal transfer function be G (s). Then the perturbed system 
can be represented in the form of an unstructured additive 
uncertainty as shown in fig. 4. W a (s) represents the weighting 
transfer function and misrepresents a set of transfer functions 
with peak magnitudes less than or equal to one for all 
frequencies. 



R(s) , 
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Fig.4 Closed Loop Uncertain System 

G(s) = Go(s) + W a (s) . A(s) 

if \A(s)\ < 1 then, \ W.(s) -A(s)\ < \ W a (s)\ 



(1) 



The weighting transfer function is chosen such that the 
following condition is satisfied. 



| W a (s) | > Max\G k (s)-G (s) \ for all co €{0, oo] 
where k = 1, 2, 3, .... n 



(2) 



where Gk(s) represents the set of n stable transfer functions 
describing the uncertain system. 

The characteristic equation of the closed loop uncertain system 
is given by, 1 + C(s) G(s) — 

Using equation (1) in the above equation we have, 



1 + C(s) \Go(s) + W a (s) . A(s)] = 

C{s)G D {s) lV Q (s).A(s) 



[1 + C(s)G D (s)] 



1 +■ 



- 



l + C(s)G D (s)" C c (s) 

The Closed Loop Nominal transfer function is given by, 
G N cl(s) = C(s) Go(s) I [1 + C(s) GofsJ] (3) 
[1+ C(s) GofsJ] [1 + G N cl(s). W a (s)A(s) /G (s)} = (4) 

Let us assume that the nominal closed loop system is stable. 
This implies that the nominal characteristic equation, 
1 + C(s) G (s) = is stable. 

The Necessary and Sufficient condition for the closed loop 

system to be stable is that equation (4) should be stable which 

implies that 1 + C(s) Go(s) — should be stable and 

[1 + G NCL (s).W a (s) A(s)/G (s) ]=0 must satisfy Small Gain 

theorem, i.e., the magnitude of the loop gain must be less than 

one. 

I G N cl(s). W a (s) A (s) I G (s) I < 1 for all co€[ 0, oo] 

Consider the worst case value of A(s) i.e., \A(s)\ — 1 for 
all co £[0, oo] in the above equation. 



I G Na (s). W a (s)/G (s) I < 1 
I Gncl(s)\ < \G (s)l W a (s)\ 



(5) 



Hence we need to make a suitable choice of Gncl(s) so that the 
above condition is satisfied. 

Consider (3). 

C(s) Go(s)l\l+ C(s) GofsJ] = G N cl(s) 
C(s)= G N cl(s) I { G (s) [1- G NCL (s)] } 

Each of the transfer functions in the above equation are 
expressed as ratios of pertinent polynomials. 



CM 



m 00 



'JfCM 



g £lCL M 



The Controller transfer function is given by, 

C C (S)/ C m (s)=G om (s) G N CLc(s)/{Goc(s)[G N CLm(s)-G NC Lc(s)]} 

(6) 

Expressing the transfer functions in equation (3) as ratios of 
pertinent polynomials, 



JNCLc I ' GtfCLm' 



- [ Cc Gqc I Cm Goin] / [ 1 ~t~ Gc Goc / Gm Goin] 
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Gnclc(s)/ G N cLm(s) = G 0c (s)/[ C m (s) G 0m (s)/ C c (s) + Goc(s)] 
Gnclc(s)/ G N cun(s) = G 0c (s)/[ C m (s) Pom(s)+G c(s]\ 
Equating the numerators and denominators on both sides, 
Gnclc(s) = Goc(s) (7) 
G NC Lm(s) = C m (s) Pom(s) + G 0c (s) (8) 
where P 0m (s) = G 0m (s) / C c (s) (9) 
Using (7) in (6) the controller transfer function is given by, 

C C (S)I C m (s) = Gom(s) /[G NC Lm(s) " G 0c (s)\ 

The general transfer function for a PID controller is given by, 

C c (s)/Cm(s)= N(s)/Ks= (ridtf+npS+nO/Ks 

We choose the numerator and the denominator of the 
controller transfer function as follows. 

From (9), G om (s) = Pom(s) C c (s) 

From (8), G NC Lm(s)-Goc(s)= P 0m (s)C m (s) 

Hence, the controller transfer function becomes, 

Gom(s)/[G NC Lm(s) -Goc(s)] = P m(s)C c (s)/ P m(s)C m (s) 

Equating this to the general transfer function for a PID 
controller, 

Go m (s)l{G N cLm(s)-Goc(s]\ = P 0m (s)N(s)/ P 0m (s).K.s 
Therefore, we have, 

Gom(s) = P 0m (s)N(s) = P Om (s)(n d s 2 +n p s+n0 (10) 

G NC Lm(s) - Goc(s) = Pom(s) .K.S (11) 

In (10), we choose the degree of Pom(s) such that it is equal to 
the degree of G om (s) minus the degree of N(s). Hence P om (s) 
and N(s) can be determined. 



From (7) and (11), 

Gnclc(s) /G NC Lm(s) = G 0c (s)/[Pom(s).K.s + Goc(s)\ 



(12) 



Now we have expressed Gncl as a function of K. That value of 
A'is chosen so as to satisfy the condition given by equation (5). 

Once A'is known, the controller parameters can be determined 
as below. 



Kd =n d /K, K p =n p /K and K, =m/K 

III. Design of Robust Controller 



(13) 



The transfer functions for the Warm air-drying chamber [xii] 
for step change in power from 0-20%, 20-40% and 40-60% of 
the range are given by, 

G(s) =10.6/ (12893 s 3 + 8058 s 2 + 226.7s + 1) for 0-20% 
G(s)= 10.5/ (2078 s 3 + 1 0346 s 2 + 230. 9s + 1) for 20-40% 
G(s)= 11.4/ (40184 s 3 + 1 0506 s 2 + 255. 6s+ 1 ) for 40-60% 
The Nominal model is given by, 

G (s)= 10.833/ (18385 s 3 + 9636.67 s 2 + 237.73s + 1) 

= Goc(s)/ Gom(s) 

Let the PID controller to be designed be given by, 

C(s) = C c (s)/ C m (s) = (n d s 2 +n p s+nj)/ Ks 

Using equation (10), the PID controller parameters are 
determined as given below. 



Gom(s) = Pom(s) (n d s 2 +n p s+nj 



(14) 



The degree of P om (s) has to be equal to one. Let 
P om (s) = pis+ 1 

Equating the coefficients of respective powers of s on both 
sides of equation (14), 

n d = 9163.743; n p =235.724; n, = l; pi= 2.0063 (15) 
Using equation (2), we choose a suitable W a . 



r 



i and Gk-Gu plots 



1 



Frequency Indirect 



Fig. 5 Frequency Response plots of \W a \, \ Gi-Go\, \ G2-Go \ and 

\G 3 -G \ 

| Go/ W a \ and Gncl for different K values are plotted. That K 
value is chosen such that the robust stability condition of 
equation (5) is satisfied. We have K=120. 

Using this K value and equation (15) in equation (13), the 
controller parameters are determined. 
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Fig. 6 Frequency Response plot of IGO / Wal and G NC l 
IV. Results 

Fig. 8 shows the closed loop step response plot with the 
introduction of the controller. It can be seen that the designed 
PID controller results in the Robust operation of the plant 
under 0-60% of the working range. 



Working 
range 


Peak value 


Overshoot 
(%) 


Rise time 
(seconds) 


Settling time 
(seconds) 


Before 


After 


Before 


After 


Before 


After 


Before 


After 


0-20% 


1.199 


1 


30.473 





36.601 


18.829 


276.473 


48.017 


20-40% 


1.213 


1.012 


32.000 


1.155 


41.752 


25.380 


323.358 


39.993 


40-60% 


1.264 


1.010 


37.535 


0.994 


37.489 


16.940 


373.362 


26.068 



Table 1 . Time response specifications before & after the 
introduction of the controller 
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Fig. 7 Closed Loop Step response without Controller 



With the introduction of the PID Controller, the maximum 
peak overshoot, maximum rise time and maximum settling 
time are reduced to 1.155%, 25.38 seconds and 48.017 seconds 
respectively. Also the steady state error is reduced from 0.917 
to which indicates an improved performance. 

V. Conclusion 

The PID controller design methodology which is simple and 
easy to use can be successfully applied to a warm-air drying 
chamber which is an example for an uncertain system. The 
results show that the designed PID controller results in Robust 
operation when the temperature is changed between and 60 
% of the maximum working range for the three operating 
points. 
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Fig. 8 Closed Loop Step response with Controller 
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